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The very first stars likely formed from metal-free, molecular hydrogen-cooled gas at the center 
of dark matter minihalos. Prior to nuclear fusion, these stars may have been supported by dark 
matter heating from annihilations in the star, in which case they could have grown to be quite 
massive before collapsing to black holes. Many remnant black holes and their surrounding dark 
matter density spikes may be part of our Milky Way halo today. Here we explore the gamma-ray 
signatures of dark matter annihilations in the dark matter spikes surrounding these black holes 
for a range of star formation scenarios, black hole masses, and dark matter annihilation modes. 
Data from the Fermi Gamma-Ray Space Telescope are used to constrain models of dark matter 
annihilation and the formation of the first stars. 
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1. Introduction 

The very first generation of stars, known as Population III.l, likely formed from the pristine 
gas at the centers of ~ lO^M© dark matter minihalos at z > 10 [1]. The response of the dark matter 
in a minihalo to the formation of a compact baryonic object at its center is a contraction of the 
dark matter density profile in and around the object. When the first stars, expected to be > lOOM©, 
ended their lives by collapsing to black holes ^, each remnant remained surrounded by a region of 
enhanced dark matter density, which we call a dark matter spike. In fact, if the first stage of stellar 
evolution is a Dark Star (DS) phase, during which the star is powered by dark matter annihilations, 
the first stars would have grown to be even larger, leaving correspondingly larger black holes and 
surrounding dark matter spikes. 

Many of these spikes, remnants of the formation of the first stars, may have merged into 
our Galactic halo, constituting Milky Way dark matter substructure today. Here we report on the 
results of Ref. [3]: the gamma-ray flux from dark matter annihilations in the dark matter spikes in 
our Galactic halo, and how we can use data from the Fermi Gamma-Ray Space Telescope (FGST) 
to constrain models of Population III. 1 star formation and/or dark matter annihilation. 

2. The First Stars and Their Dark Matter Spikes 

Population III. 1 stars could only form when a cooling mechanism for the collapse of the bary- 
onic cloud arose. The first accessible mechanism to cool the gas was via excitations of molecular 
hydrogen, the fraction of which present in a minihalo is related to the temperature, which in turn 
can be written in terms of the mass and redshift of the minihalo. We use the parametrization of 
Ref. [4] for the minimum halo mass in which star formation could occur: 



We take the maximum halo mass for Population in.l star formation to be = XQPMq, though 
our results are not sensitive to this choice. 

At some time between the beginning of Population III. I star formation and the end of reion- 
ization at z 6, massive Population III.l (and DS) formation must have given way to subseqent 
formation of less-massive stars [5], however there are few constraints on when this transition oc- 
cured. Here we consider three scenarios for the termination of Population III. 1 star formation at 
redshift Zf. These scenarios are hereafter noted as Early, Intermediate, and Late, following Ref. [5], 
withz/Ri23, 15, and II, respectively. For each case, we assume that Population HI. I star formation 
was possible in any minihalo with a mass between M^"^" and at redshift z>Zf- 

In fact, not every minihalo meeting the above criteria must have hosted a Population III. 1 star. 
We therefore parametrize the fraction that did as fjjs- Neglecting black hole mergers, a discussion 
of which can be found in [3], the comoving number density of black holes as a function of redshift 
is then 



Stars in the mass range ~ 140 — 260Mq would have ended their lives as pair instabiUty supemovae, leaving no 
remnants [2]. We do not consider these objects here. 




(2.1) 



Nbh{z) = fosNhaloiz), 



(2.2) 



2 



Signatures of Dark Star Remnants in the Galactic Halo 



Pearl Sandick 



where Nuaio (z) is the comoving number density of minihalos in which Population III. 1 star forma- 
tion was possible. 

Assuming that each viable minihalo hosted a Population III.l star (Jds = 1), we determine the 
z = distribution of dark matter spikes throughout the Galactic halo from the Via Lactea II (VL- 
II) cosmological N-body simulation [6]. In Fig. [l| we show the number densities of dark matter 
spikes inside the Milky Way halo as functions of Galactic radius for Early, Intermediate, and Late 
Zf. In the Early case, Pop. III.l star formation terminates at z 23, so there were the fewest stars, 
and therefore the fewest black holes and surviving density spikes today; 409 spikes in our Galactic 
halo. In the Intermediate and Late cases, we find 7983 and 12416 dark matter spikes in our Galactic 
halo, respectively. In a similar analysis, Ref. [7] found the distribution of spikes given by the grey 
curve in Fig. [l], resulting in 1027 spikes in our Galactic halo. For comparison, the total dark 
matter density profile at z = in VL-II is also shown; although the normalization of these points 
is arbitrary, it is useful to illustrate that the total dark matter profile is more extended than the 
distribution of black holes with dark matter spikes. 
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Figure 1: The number density of black hole spikes in the Milky Way as a function of Galactic radius 
for star formation models with Early (green). Intermediate (red), and Late (blue) z/ as described in 
the text. The black points illustrate the total dark matter density profile at z = 0. Also shown as a 
solid grey curve is the analytical fit found in Ref. [7, 8]. 

The density profile of an individual dark matter spike surrounding a black hole is determined 
by adiabatic contraction of the dark matter halo around the central mass. We assume Navarro, 
Frenk, and White (NEW) profiles for both the baryons and dark matter, and compute the contracted 
dark matter profile using the Blumenthal et al. prescription for adiabatic contraction [9]. The 
resulting halo profiles are roughly power-law in nature with a high density cut-off imposed to 
account for the annihilation of the dark matter in the spike; p^ax = "^;t:/((<7v)?g//), where nij^ is the 
mass of the dark matter particle, (av) its annihilation cross section times velocity, and Ibh is the 
lifetime of the central mass, roughly 1.3 x 10^° years for a star that formed at z = 15. 

While standard Population III. 1 stars are expected to have masses of ~ IGGM©, if dark matter is 
capable of self-annihilating, then the first stars may have been powered by dark matter annihilations 
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for some period of time prior to nuclear fusion. This first phase of stellar evolution is known as the 
Dark Star (DS) phase [10] and may have lasted anywhere from a few hundred thousand years to 
millions, or even billions, of years. During the DS phase, the star remains cool enough to continue 
to accrete baryonic matter, and may grow to ~ lOOOM© [11] or even as large as > IQr'Mr.y [12], 
depending on the details of how the dark matter is depleted and replenished in the star. 

When the dark matter fuel inside a DS runs out, it collapses and heats up to become a standard 
fusion-powered star, which, at the end of its life, will likely undergo core collapse leaving a black 
hole remnant. We use the potential existence of the DS phase to motivate consideration of black 
holes as large as IO^Mq. 



3. Gamma Ray Signal from Dark Matter Annihilations 

For a Majorana dark matter particle with mass mj(^ and annihilation cross section times velocity 
(av), the rate of WIMP annihilations in a dark matter spike is 

r=^/ drAnr^plMir), (3.1) 



nun 



with r,„in and rmax defining the volume of the dark matter spike in which annihilations occur. 

We choose as a benchmark scenario (av) = 3 x 10^^^ cm^s^\ in agreement with the mea- 
sured dark matter abundance today for thermal WIMP dark matter, and consider several WIMP 
candidates, defined by mass and annihilation channel. We consider annihilations of WIMPs with 
masses from 100 GeV to 2 TeV to Standard Model final states bb, W+W^, T+T~, and }X^}X^. The 
resulting spectrum of photons from annihilation to final state / is dNf/dE. For %% > the 

photon spectrum comes only from final state radiation. In the following analysis, we assume that 
the branching fraction to each final state is 1 , though, in principle, some combination of final states 
is possible. 

The differential flux of neutral particles from annihilations to final state / in a dark matter 
spike with radius r^ax located some distance D from our Solar System is given by 

d^^J^dN^ 
dE 47iD^ dE' 

forD > 

^max- If ^ single spike is a sufficiently bright and compact source of gamma-rays, it may 
have been identified as a point source and recorded in the FGST First Source Catalog [13]. Spikes 
bright enough to be identified as point sources must not be brighter than the brightest source in the 
the FGST catalog. Requiring that the flux not exceed this brightness establishes a minimal distance, 
^mfn' beyond which the spike must be located. Similarly, we can define a maximal distance, D^^^, 
as the distance beyond which a spike would not be bright enough to have been detected at 5 a 
significance (see [3] for details). All spikes which could have been identified as point sources lie 
between DfJ, and D^,. 

In Fig. ^ we display in the left panel the minimal distance, D^f„, at which a point source may 
be located in order not to exceed the largest flux from any point source measured by FGST for eight 
example WIMP annihilation cases. For most of the dark matter models shown, if the central black 
hole is < 10^ M0, spikes may be located within 1 kpc of our Solar System. However, for all choices 
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Figure 2: In the left panel, we display the minimal distance from our Solar System of a single spike 
such that it does not exceed the flux of the brightest source in the FGST First Source Catalog, 
as a function of central black hole mass for Intermediate Zf- From top to bottom, the contours are 
for 100 GeV WIMPs annihilating to bb (blue dashed), 100 GeV WIMPs annihilating to W+W- 
(orange dashed), 100 GeV WIMPs annihilating to t+t^ (magenta dashed), 1 TeV WIMPs annihi- 
lating to bb (blue sohd), 1 TeV WIMPs annihilating to W+W^ (orange sohd), 100 GeV WIMPs 
annihilating to (green dashed), 1 TeV WIMPs annihilating to T+T^ (magenta), and 1 TeV 

WIMPs annihilating to (green solid). In the right panel, we show the maximal distance from 

our Solar System of a single spike such that it would appear as a > 5a point source to Fermi, D^^, 
for the same cases. The horizontal black line indicates the distance to the Galactic center. 



of Zf considered here, we expect less than one dark matter spike within 1 kpc of our Solar System, 
even for f^s = 1. In the right panel of Fig. ^ we show the maximal distance, D^^^, at which a 
single dark matter spike would have been identified by FGST as a > 5a point source. For 100 GeV 
WIMPs annihilating to bb or IV+W^ around large black holes, as in the upper right portion of the 
plot, all dark matter spikes in the Milky Way halo would have been identified as point sources by 
FGST. 

Many of the dark matter spikes in the Milky Way halo may be too faint to be identified as 
point sources by FGST and would therefore contribute to the diffuse flux. Based on the simulated 
distribution of Population III.l/DS remnant spikes in the Milky Way halo from VL-II, we calculate 
the contribution to the diffuse gamma ray flux from all spikes that are decidedly too faint to have 
been identified as point sources by FGST. These results are presented in Fig. ^ for annihilations of 
100 GeV WIMPs to bb (left panel) and (right panel) along with the FGST-measured diffuse 

gamma-ray flux [14]. We see that the all-sky averaged diffuse flux is well below that measured by 
FGST in all models shown at low energies, but can rise to the level of the FGST-measured flux at 
intermediate and higher energies. 

4. Constraining //)5 

To this point we have focused only on /ds = 1 , though in fact /ds may be ^ 1 . It is possible 
to constrain fi,s in two ways: First, from the FGST measurement of the diffuse gamma-ray flux, by 
requiring that the diffuse flux from dark matter annihilations around spikes in the Milky Way halo 
not exceed the measured flux in any of the nine FGST energy bins in Ref. [14] by more than 3a. 
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Figure 3: Diffuse gamma-ray flux from dark matter spikes for annihilations of 100 GeV WIMPs to 
bb (left) and pi^pL^ (right). In each panel, Early (green, dotted). Intermediate (red, solid), and Late 
(blue, dashed) z/ are shown for central black holes of mass lO^M© (thick curves) and 10^ M0 (thin 
curves). The black points represent the FGST-measured diffuse gamma-ray flux [14]. 



The resulting maximal values of /05 are presented in Fig. Q for annihilations of 100 GeV WIMPs 
to bb (left panel) and (right panel) as the open points, with Early, Intermediate, and Late Zf 

represented in green, red, and blue, respectively. 

A second way to constrain fos is by requiring that one should expect to find less than one dark 
matter spike within D^,f„, the minimal distance at which a spike may be located such that it is not 
brighter than the brightest point source in the FGST First Source Catalog [13]. The maximum f^s 
as found in this method is displayed in Fig. ^ as the solid points in each scenario. These are the 
most conservative limits on f^s, as they are based on the flux not exceeding that of the brightest 
FGST point source, known to be the Vela pulsar. If we require that the flux from a dark matter 
spike not exceed that of the brightest unassociated point source, stronger limits on f^s would be 
obtained, though there is no gaurantee that there is not a dark matter spike along our line of sight 
to Vela. 
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Figure 4: Maximum fj)s as a function of central black hole mass for 100 GeV WIMPs annihilating 
to bb (left) and jJ-^jJ.^ (right). Green circles, red squares, and blue diamonds are for Early, Inter- 
mediate, and Late z.f, respectively. The solid markers are the limits from point source brightness, 
while the open markers are from the diffuse gamma-ray flux, as described in the text. 
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5. Conclusions 

We have show how FGST gamma-ray observations can be used to place limits on the properties 
of dark matter spikes around black holes in our Galactic halo. We compare the gamma-ray flux 
from dark matter annihilations around spikes to both the gamma-ray flux from point sources and 
the diffuse flux observed by FGST, and constrain the maximum fraction of viable minihalos that 
could have hosted Population III. 1/DS star formation. In general, it is clear that the bounds are the 
strongest for the largest black hole masses and if star formation persisted to low redshift (in which 
case the total possible number of dark matter spikes in our Galactic halo is largest). These limits 
also depend sensitively on the dark matter annihilation channel, as one can see by comparing the 
results for annihilations to bb and ix^jx~ final states. If/when the annihilation properties of particle 
dark matter become known, we may be provided with some interesting hints about the formation 
of the first stars. 
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